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ABSTRACT 



Context. The surface layers of the Sun are strongly stratified. In the presence of turbulence with a weak mean magnetic 
field, a large-scale instability resulting in the formation of non-uniform magnetic structures, can be excited over the 
scale of many turbulent eddies or convection cells. This instability is caused by a negative contribution of turbulence 
to the effective (mean-field) magnetic pressure and has previously been discussed in connection with the formation of 
active regions and perhaps sunspots. 

Aims. We want to understand the effects of rotation on this instability in both two and three dimensions. 
Methods. We use mean-field magnetohydrodynamics in a parameter regime in which the properties of the negative 
effective magnetic pressure instability have previously been found to be in agreement with those of direct numerical 
simulations. 

Results. We find that the instability is suppressed already for relatively slow rotation with Coriolis numbers (i.e. inverse 
Rossby numbers) around 0.2. The suppression is strongest at the equator. In the nonlinear regime, we find traveling 
wave solutions with propagation in the prograde direction at the equator with additional poleward migration away from 
the equator. 

Conclusions. The prograde rotation of the magnetic pattern near the equator is argued to be a possible explanation for 
the faster rotation speed of magnetic tracers found on the Sun. In the bulk of the domain, kinetic and current helicities 
are negative in the northern hemisphere and positive in the southern. 

Key words, magnetohydrodynamics (MHD) - hydrodynamics - turbulence - Sun: dynamo 
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1. Introduction 

In the outer parts of the Sun, energy is transported through 
turbulent convection. The thermodynamic aspects of this 
process are we ll understood through mixing length theory 
(|Vitensel . ll953D . Also reasonably well understood is the par- 
tial conversion of kinet i c energy into magnet i c ener gy via 
dynamo action (I ParkeiiUQTOHZeldovich et al.l ll983'). Most 
remarkable is the possibility of generating magnetic fields 
on scales much larger and much longer than those of the 
turbulence itself. This has now been seen in many three- 
dimensional turbulence simulations (jBrandenburd . l200lt 
[Brandenburg fc Subramanianl . |2005[ ). but the physics of 
this is best understood in terms of mean-field theory, which 
encapsulates the effects of complex motions in terms of ef- 
fective equations for mean flow and mean magnetic field 
([Mofratt, 1978; Parker, 1979; Krausc k Radlcr, 1980). 

The effects of stratification are usually only included 
to leading order and often only in connection with rota- 
tion, because the two together give rise to the famous a 
effect, which is able to explain the generatio n of large- 
scale magnetic fields (jKrause fc Radled . 119801) . In recent 
years, however, a completely different effect arising from 
strong stratification alone has received attention: the sup- 
pression of turbulent pressure by a weak mean mag- 
netic field. This effect mimics a negative effective (mean- 
field) magnetic pressure due to a negative contribution of 



turbulence to the mean magnetic pressure. Under suit- 
able conditions, this leads to the negative effective mag- 
netic pressure instability (NEMPI) which can cause the 
formation of magnetic fiux concentrations. In turbulence 
simulations, this in s tabilit y has been seen only recently 
(jBrandenburg et "all . 1201 ll ). because significant scale sep- 
arati on is needed to overcome t he effects of turbulent diffu- 
sion (jBrandenburg et"aI1 . l2012h . Mean- field considerations, 
howe ver, have predicted the existence of NEMPI for a long 
time (iKleeorin et al 



1989[ | 1990HKleeorin fc Rogachevskii , 



'1994': 'Kleeorin et al.'.'l996':'Rogachevs kii fc Kleeorinl . 12007 
Brandenburg, Klccorin, & Rogachcvskij , 2010t) 



Despite the potential importance of NEMPI, many ad- 
ditional effects have not yet been explored. The idea is that 
NEMPI would interact with the global dynamo producing 
the large-scale magnetic field for NEMPI to act upon. Thus, 
the field needs to be self-consistently generated. Ideally, 
global geometry is needed, and such calculations should be 
three-dimensional (3-D), because one expects flux concen- 
trations not to be two-dimensional (2-D) or axisymmetric. 
New mean-field coefficients will appear in such a more gen- 
eral case, and not much is known about those coefficients. 
Nevertheless, although other terms may appear, it will be 
interesting to investigate the evolution of NEMPI in more 
realistic cases with just the leading term responsible for the 
instability. 
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The purpose of the present paper is to include the ef- 
fects of rotation as it appears in a local Cartesian domain 
at a given latitude in the Sun. Our main goal is to deter- 
mine the dependence of growth rate and saturation level of 
NEMPI on rotation rate and latitude, and to characterize 
rotational effects on the resulting flux concentrations. For 
the purpose of this work, we restrict ourselves in the follow- 
ing to a mean-field treatment and denote averaged quanti- 
ties by an overbar. Furthermore, we make the assumption 
of an isothermal equation of state. This is of course quite 
unrealistic, as far as applications to the Sun is concerned. 
However, it has been found earlier that NEMPI has similar 
properties both for an isothermal layer with an isothermal 
equation of state as well as a n early isentropic one w ith the 
more general perfect gas law (jKapvla et al.l . |2012[ ). Given 
that our knowledge of NEMPI is still rather limited, it is 
useful to consider the new effects of rotation within the 
framework of the conceptually simpler case of an isother- 
mal layer. 

One of the remarkable insights is that NEMPI can occur 
at any depth, depending just on the value of the mean mag- 
netic field strength. However, for a domain of given depth 
extent, the instability can only occur if the location where 
the dependence of effective turbulent pressure on the ratio 
field strength to equipartition value has a negative slope. 
Once this is obeyed, the only other necessary condition for 
NEMPI to occur is that the turbulent diffusivity is small 
enough. In practice this means that there are sufficiently 
many turbulent edd i es wi t hin the domain of in vestigation 
(|Brandenburg et all . 120121: iKemel et al.l l2012b[ ). 

We begin with the model equations, present then first 
the linear theory of NEMPI in the presence of rotation, 
and consider then two- and three-dimensional numerical 
models. 



on the right hand side of Equation ^ determines the tur- 
bulent contr ibution to th e mea n Lorentz for ce. Following 
Brandenburg et al.l (|201^ and i Kemel et all (pbl2a) . the 
function qp{/3) is approximated by: 

,p(« - (6) 

where l3^, and /3p are constants, /? = B/Bcq is the mod- 
ulus of the normalized mean magnetic field, and B^q — 
y^JtopUrms is the equipartition field strength. The angular 
velocity vector Jl is quantified by its scalar amplitude fl 
and colatitude 9, such that 

fl = n{-am9,0,cose) . (7) 

In this arrangement, z corresponds to radius, x to colati- 
tude, and y to azimuth. 

Following the simplifying assumption of recent di- 
rect n umerical simulations of NEMPI (Brandenburg et "all 
l201lD . we assume that the root-mean-square turbulent ve- 
locity, Mrms, is coustaut in space and time. For an isothermal 
density stratification, 

p = PQexp{-z/Hp), (8) 

where Hp = c^/g is the density scale height, we then have 
Bcq{z). To quantify the strength of the imposed field, we 
also define Boqo = Bcq{z = 0). The value of Wrms is also 
related to the values of rjt and i>t, which we assume to be 
equal, with rjt = I't = Wr,ns/3fcf , where fcf is the wavenumber 
of the energy-carrying eddies of the underlying turbulence. 
This formula assumes that the relevant correlat ion time is 
(wrmsfcf , which turns out to be fairly accurate ()Sur et al.l 
,2008.) . 



2. The model 



3. Linear theory of NEMPI with rotation 



We consider here an isothermal equation of state with con- 
stant sound speed Cg , so the mean pressure isp = pc^ . The 
evolution equations for mean velocity U, mean density p, 
and mean vector potential A, are 



^ = -2^ xU- c^V Inp + fif + + ^K, 



Dp 

m 

dA 
'dt 



= UxB-{r]t+ r?) J, 



(1) 
(2) 

(3) 



where D/D< = d/dt + U ■ V is the advective derivative, 
g = (0, 0, —g) is the acceleration due to the gravity field. 



= [vt + v) (V2J7 + i VV • U + 2SV Inp) 



(4) 



is the total (turbulent plus microscopic) viscous force, and 
Sij = ^(/7i,j-l-(7j,i) — ^i5ijV -C/ is the traceless rate of strain 
tensor of the mean flow. The mean Lorentz force, is 
given by 



(5) 



where J — x B/po is the mean current density, po is 
the vacuum permeability, and the last term, iV((7pB^), 



In this section we study the effect of rotation on the growth 
rate of NEMPI. Follo wing earlier work (e.g., the appendix of 
IKemel et all l2012bD . we neglect for simplicity dissipation 
processes, use the anelastic approximation, V • 'pU = 0, 
and assume that the density scale height Hp = const. We 
consider the equation of motion, ignoring the U ■'VU non- 
linearity. 



dU{t, X, z) 
di 



-20 X U 



= Vptot 
P 



g, 



(9) 



where ptot = p + PcS is the total pressure consisting of the 
sum of the mean gas pressure p, and the effective magnetic 

pressure, PeS — {1 — qp)B /2, where B — \B\. Here and 
elsewhere the vacuum permeability is set to unity. Let us 
assume for simplicity that dy = 0, and that the mean mag- 
netic field has only a y-component, B ~ (0, By{x, z),0), so 
the mean magnetic tension, B ■ VB in Equation Q van- 
ishes. 

Taking twice the curl of Equation ([5]), and noting fur- 
ther that i: • V X V X C7 = -AUz + V • U, we obtain 



V,{U -V In p)] = ■ V(V x U), 



Ptot \ ^xP 

p 



- V 



Ptot 

p 



"^zP 



(10) 
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where we have used the anelastic approximation in the form 
V-C/ = — J7-Vlnp and the fact that under the curl the 
gradient can be moved to p. We have a lso taken into account 
that ny = and have used Eq. (30) of iKemel et"all (|2012b[ ) 
to relate the double-curl of (Vptot)/p to the last term in 
Equation ([TU]). The first term on the right hand side of 
Equation dTU]) for Uz is proportional to (V x U)z- Taking 
the z component of the curl of Equation ^ we obtain the 
following equation for (V x U)z- 



|(Vx[/).^2(o.V-§^)c/.. 

The induction equation for By{x, z) is given by 



(11) 



(12) 



where D/Dt = d/d + (7 • V is the advective derivative. 
Note that for a magnetic field with only a y-component, 
but d/dy — 0, there is no stretching term, so there is no 
term of the form B ■ Vt/. 

We linearize Equations p^ -([T ^ . indicating small 
changes by 5. We consider an equilibrium with a constant 
magnetic field of the form (0, i?o,0), a zero mean velocity 
and the fluid density is given by Equation dS]). We take into 
account that the function q-p = q-piP) depends bo th on B 
and on p, which implies that (jKemel et al.l . l2012bf ) 



Ptot 



1 2 M dgp 



while 



Ptot 



hv\[l-qp- 



dln/32 



JBy Sp 



P 



1 



P J ' -' \ " dln/32 J Hp' 
The linearized system of equations reads: 



(13) 



(14) 



dt 



Hn 



Hp d/32 So 
-20 • V(V X 6U)z 



(15) 



|:(V X 5U)z =2(n-V-^\ 6Uz, (16) 



dt 



dSBy 

at 



-Bn 



SUz 



(17) 



where Vcs{f3) = 5 [1 ~ 9p(/^)] /^^ the effective magnetic 
pressure normalized by the local value of i?oq- 

Introducing a new variable Vz — y/pSUz in 
Equations P5|) - PT)) and after simple transformations we 
arrive at the following equation for one variable Vz'- 

(18) 



where 

Xliz) = -2 



vl{z) dPeff(^) 

HI d/32 ■ 



In the WKB approximation, which is valid when 
kz Hp ^ 1, i.e., when the characteristic scale of the spatial 



variation of the perturbations of the magnetic and velocity 
fields are much smaller than the density height length. Hp, 
the growth rate of the large-scale instability is given by 



A = 



^ 2 '^a; _ 2 
^2 '^incrt 



1/2 



(20) 



where Winoit = 217 • k is the frequency of the inertial waves. 
Here, fe = fc/fc is the unit vector of k. Note that A is 
either real or purely imaginary, so no complex eigenval- 
ues are possible, as would be required for growing oscil- 
latory solutions. Without rotation t he growth rate of the 
large-scale instability (N EMP P is (iKleeoriii et al.l . 119931: 
iRogachevskii fc Kleeoriil l2007t iKemel et al.l . l2012cl ) 



d/32 



1/2 



(21) 



(22) 



and a necessary condition for the instability is 

^<0. 

d/32 

NEMPI can be excited even in a uniform mean magnetic 
field, and the source of free energy of the instability is pro- 
vided by the small-scale turbulence. In contrast, the free 
energy in Parker's magnetic buoyanc y instabilitv (jParker , 
19 6a) or in t he interchange instability (jTserkovnikovl 1960 : 
|PriestLll982|) is drawn from the gravitational field. Both in- 
stabilities are excited in a plasma when the characteristic 
scale of variations of the original horizontal magnetic field 
is smaller than the density scale height. 

The rotation reduces the growth rate of NEMPI, that 
can be excited when k^/k > Winert/Ao and d'Pcff/d/32 < 0. 
In opposite case, k^/k < Wincrt/Ao, the large-scale instabil- 
ity is not excited, while the frequency of the inertial waves 
is reduced due to the effective negative magnetic pressure. 

For an arbitrary vertical inhomogeneity of the density, 
we seek a solution of Equation (ITS)) in the form Vz {t, x, z) = 
V{z) exp(Ai -|- ikx x) and obtain an eigenvalue problem 



where 
A2 = 



2 ' ^"^"^ tkrV.-A^kl 



1 



A2 + 417? 



,_x^-xliz) + 4nl 



AH^p 



A2 -I- 4f]2 



V{z) = 0, (23) 



(24) 



and A is the eigenvalue. Equation (^5]) can be reduced to 
the Schrodinger type equation, ^" — U{R) ^I' = 0, via the 
transformation 



*(i?) = VRV{z) exp (i^^^k^z 



R{z) 



AO 



u2 /32 
rmsr^p 



(25) 
(26) 



where vaq = B^l^fPo is the Alfven speed based on the 
averaged density, the potential U{R) is 



(19) U{R) 



rA2 /A2 + 4f}^ 



i?(A2 + 4f^2) 



R \y + mi 



(27) 
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0.02, 61 = and 17 = 0.01, 0.1, 



and we have used Equation ^ for qp with jS^, — /3p y/Qpo and 
QpO — QpiP — 0). As follows from Equation (l?7l) . the poten- 
tial, U (R), is positive for i? — and R — oo. Therefore, for 
the existence of the instability, the potential should have a 
negative minimum. This is possible, when qpQ > {1 + R)^. 
When the potential U{R) has a negative minimum, there 
are two points i?i and R2 (the so-called turning points) 
in which U{R = R1.2) — 0. Figure [T] shows tanh t/(i?) for 
different values of f2. This representation allows us to dis- 
tinguish better the behavior for small values of U (R) . 

Using Equation (P7|) and the condition U[R — i?i,2) — 
0, we estimate the maximum growth rate of the instability 
as 



V2 



A^ - + - An^f + mlxl] 



1/2 



1/2 



where 



A* 



/?,W™s [i?li?2(2 + i?l +i?2)] 



1/2 



{I + Ri){l + R.2) 



(28) 



(29) 



By defining a = AVi^/Xl, Equation ^ can also be written 
as 



1 



A/A, = ^ 1 - (7 + (1 - 2(7 sin^ 61 -I- CT^) 
V2 L 



1/2 



(30) 



Note that for cr ^ 1, we obtain that A/ A* = cos6', which 
is independent of the value of a. In Figure [2] we plot the 
dependence of A/ A, on 6 for different values of a. This 
analytic result shows that, if A, were independent of 9, the 
growth rate remains unaffected by rotation at the poles, 
and is suppressed near the equator for VL — 0(A,). 

Unfortunately, the analytic result does not allow an easy 
inspection concerning the sign of the frequency and hence 
the direction of the travelling wave solutions. Therefore we 
turn in the following to numerical simulations of the full 
2-D and 3-D mean-field equations. 

4. Numerical Results 

In this section we discuss numerical mean-field modelling. 
We consider a computational domains of size or Li^ 
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Fig. 2. Theoretical dependence of A/A, on for different 
values of a. 



with periodic boundary conditions in the horizontal direc- 
tion(s) and stress-free perfect conductor boundary condi- 
tions in the vertical direction. The smallest wavenumber 
that fits horizontally into the domain has the wavenum- 
ber fci = 27r/L. The numerical simulations are performed 
with the Pencil CodeQ, which uses sixth-order explicit 
finite differences i n space and a third-order acc urate time 
stepping method (jBrandenburg fc Dobled . [20021 ) . As units 
of length we use k^^ and time is measured in units of 
(cgfci)"^. For the models presented below, we use ^po = 20 
and Pp = 0.167, which corresponds to /3* = 0.75, and is 
appropriate for the parameter regime in which i?m ~ 18 
and fcf/fci = 30 (Kcmcl ct al., 2012c). 

An important non-dimensional parameter is the Coriolis 
number, Co = 2f2/ui-insfct. Using fcf = Wrms/377t, we can 
express this in terms of the parameter Cn — ^/rjtkf, which 
is often used in mean-field dynamo theory. Thus, we have 



Co = 67;tf7/w?ms = 6 {ritki/u^msTCn- 



(31) 



Motivated by the analytic results of the previous section we 
normalize the growth rate of the instability alternatively 
by a quantity A,o = P*u^ms/ Hp. In the following we take 
Wims/cs = 0.1. Further, we use vt — fit — 10~^Cs/fcf, so 
that kiHp « 33 and ?7tfci/urms = 10^^. This also means 
that for Cl — 0.01, for example, we have 2i7/A,o — 0.27 and 
Co = 0.006. 

We use either Bo/Bcqo ~ 0.1 or 0.05. We recall, how- 
ever, that the growth rate does not depend on this choice, 
provided the bulk of the eigenfunction fits into the domain, 
which is here the case for both values of Bq . (For the smaller 
value of Bq the maximum of the magnetic structures (i.e., 
the maximum of the eigenfunction in z) is slightly higher up 
in the domain, but in both cases the maximum is contained 
within the domain.) 

We discuss first the and 9 dependence of 2-D and 3-D 
solutions. Using 6* = 0°, 45°, and 90°, corresponding to 90°, 
45°, and 0° latitude, we find that NEMPI is suppressed for 
rotation rates around ft ~ O.OlCgfci and 0.025 in 2-D and 
3-D, as can be seen in Figures [3] and |31 This corresponds 
to Co = 0.006 and 0.015, which are remarkably small val- 
ues. Note the similar behavior in 2-D and 3-D: NEMPI is 
suppressed for even smaller values of 2^l/X^,o as 9 increases. 



^ |http: //pencil-code .googlecode . coiii| 
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Next, we vary 9. As expected from the results of 
Section [21 and as already seen in Figures [3] and SI the largest 
growth rates occur at the poles {9 = 0), and NEMPI is most 
strongly suppressed at the equator. The growth rate as a 
function of 6 is given in Figure [S] for two values of 2ri/A*o, 
showing a minimum at 9 = 90° (i.e., at the equator). In 
the upper panel of Figure [5l we have used two-dimensional 
results, i.e. we have restricted ourselves to solutions with 
d/dy = 0, as was also done in Section [31 However, this is 
only an approximation to the fully three-dimensional case. 
The usefulness of this restriction can be assessed by com- 
paring two-dimensional and three-dimensional results; see 
the lower panel of Figure \5\ While the 9 dependence is 
roughly similar in the two- and three-dimensional cases, 
the growth rates are by at least a factor of two smaller in 
the two-dimensional case. 

To determine the oscillatory frequency of the growing 
mode, we consider the values of Uy{xi,t) and By{xi,t) at 
a fixed point To determine the oscillatory frequency, we 
consider the values of Uy{xi,t) and By{xi,t) at a fixed 
point Xi within the domain. As can be seen in Figures [SI 



Fig. 5. Dependence of A/ A*o on 6* for two values of 2r2o/A*o 
in 2-D (upper panel) and comparison of 2-D and 3-D cases 
(lower panel). 

and [71 their frequency and amplitude depend both on il 
and 9. The oscillations are not always harmonic ones, and 
can be irregular with variable periods, making the period 
determination more difficult. Nevertheless, the frequencies 
for Uy and By are similar over broad parameter ranges, 
as is expected from linear theory. For f2o/A*o > 0.25 at 
9 = 60°, NEMPI is no longer excited, but there are still 
oscillations in Uy{xi,t) which must then have some other 
reason. 

We find a substantial variation of the amplitude for the 
values of maximum growing rate for ^l = 0.01 and = 0.02. 
(Note that the high frequency in Uy and By in Figure [HI 
corresponds to a random small-amplitude change.) The fre- 
quency of the oscillations is very small at the poles, but it 
reaches a maximum at = 45 and decreases again toward 
the equator. 

Given the combined presence of rotation and stratifica- 
tion, we expect the resulting velocity and magnetic fields 
to be helical. We plot in the upper panel of Figure [51 rel- 
ative kinetic, current, and cross helicities. These are here 
abbreviated in terms of the function 



nip,q) = {p-q)/VW)W) 



(32) 



where p and q are two arbitrary vectors. Here, (•) de- 
note xy_ averaging. Note that the relative kinetic helicity, 
H{W, U), where W = V x 17 is the mean vorticity, varies 
between nearly -1-1 in the lower part and —1 in the upper 
part. This change of sign is familiar from laminar convection 
where upwellings expand to produce negative helicity in 
the upper parts, and downwelli ngs also expand as they hi t 
the bottom of the domain (e.g. [Brandenburg et all . Il990[ ). 
However, in the lower part of the domain both U and W 
are relatively small, as can be seen by considering their 
relative amplitudes, A{U) and y^(W), where 



Aip) = {p')/{{p')), 

with ((•)) being defined as volume averages. 



(33) 
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Fig. 8. Dependence of various relative helicities and rela- 
tive amplitudes on z for the case with 9 = 0° and Co = 0.03. 
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Finally, let us turn to_thc spatial structure of NEMPI. 
In Figure [9] we compare By at different times and latitudes 
for the 2-D runs. Note that in the exponentially growing 
phase of NEMPI, the structures do not propagate (or move 
only very slowly). Traveling wave solutions occur mainly in 
later stage of NEMPI, i.e., in the saturated regime. Next, we 
consider the 3-D case. In Figure [TU] we show visualizations 
of the magnetic field on the periphery of the computational 
domain for 4 different times for 9 = 0. Note that magnetic 
structures are inclined in the xy plane. This is a direct re- 
sult of rotation. As expected, the inclination is for negative 
values of J7; see Figure [TT] The inclination angle is about 
30°, corresponding to 0.5 radians, which is not compatible 
with the value of Co « 0.03, but it is closer to the value of 
f2/A*o ~ 0.65. However, in this connection we should stress 
that we have imposed periodic boundary conditions in the 
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Fig. 9. Evolution of By in the xz plane for JIq = 0.01 
(corresponding to Co = 0.006) and Bo/Bcqp — 0.1 for 9 = 
0°, 9 = 45° , and 6 = 90° near the time when the instability 
saturates. The direction of fl is indicated in the last row. 



y direction, which means that the inclination angles change 
only in discrete steps. 

Returning to the case of positive values of fi, but 9^0, 
we note a slow migration of the magnetic pattern to the left 
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t/T^ = 2.50. 



1 1 n ' 'I ^^^B 





'i,/-S=, t/T^ = 3.50 





Fig. 10. Visualization of By on the periphery of the computational domain for 4 times (normalized in terms of T^) 
during the nonlinear stage of the instability for 9^0° (corresponding to the equator) and Co = 0.03, corresponding to 
2n/X^Q « 1.3. 
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Fig. 11. Same as Figure[TUJ but for a negative value of f2, i.e., Co ~ —0.03, corresponding to 2f2/A*o ~ —1.3. 
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Fig. 12. Visualization of By on the periphery of the computational domain for 4 times (normalized in terms of T^) 
during the nonlinear stage of the instability for 6 = 45° and Co = 0.03, corresponding to 2f2/A*o ~ 1-3. 
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Fig. 13. Visualization of By on the periphery of the computational domain for 4 times (normalized in terms of T^) 
during the nonlinear stage of the instability for 9 — 90° (corresponding to the equator) and Co = 0.013, corresponding 
to 20/A,o ~ 0.5. 



(here for 9 = 45°), corresponding to poleward migration; plane. Finally, for 9 = 90° we see that the pattern speed 
see Figure [T^l Note also that the field is tilted in the xy corresponds to prograde motion; see Figure [T31 
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5. Conclusions 

Although the physical reality of NEMPI has recently been 
confirmed by direct numerical simulations, its potential role 
for producing large-scale magnetic structures in the Sun is 
yet unclear and we are only now beginning to investigate 
its properties under conditions that are astrophysically im- 
portant. Rotation is ubiquitous and clearly important in 
the Sun. The present work has now shown that already 
for rather slow rotation the instability is suppressed. This 
is rather surprising, because rotational effects normally be- 
come important only when n is comparable with the inverse 
turnover time, which is here defined as Mrms^f- The reason 
for this behavior is probably related to the fact that the 
growth rate of the instability does also not scale with the 
inverse turnover time, but that it is closer to the turbulent 
diffusive time, which is slower than the turnov er time by 
the sq uare of the scale separation ratio ( Branden burg et ahl . 
1201 It) . However, the present work now suggests that this is 
not quite right either and that the correct answer might 
be something in between. Indeed, we find here that growth 
rate and critical rotation rate are close to the parameter 
A*o = f3i,u^nis/ Hp, which can be smaller than the afore- 
mentioned turnover time by a fac tor of 40, although i n so- 
lar convect ion, where kiHp ^ 2 .4 (|Kemel et al.l . l2012c[ ) and 

« 0.23 (|Kemel et all . l2012bD . it is estimated to be only 
« 10 times smaller. 

The suppression is strongest at the equator, where f2 
is perpendicular to the direction of the gravity field, i.e., 
fl ■ g = 0, and less strong at the poles where fl and g are 
either parallel (south pole) or antiparallel (north pole). In 
the absence of rotation, the mean magnetic field only varies 
in a plane that is normal to the direction of the imposed 
mean magnetic field, i.e., k ■ Bq = 0, where k stands for 
the wave vector of the resulting fiow and magnetic field. 
However, in the presence of rotation the orientation of this 
plane changes such that now k ■ (Bq + A~o^r2 x Bq) = 0. 

At intermediate latitudes, i.e., with angles between fl 
and g in the range from 0° to 90° colatitude, the magnetic 
field pattern propagates slowly in the negative x direction, 
corresponding to poleward migration. The significance of 
this result is unclear. Had it been equatorward migration, 
one might have been tempted to associate this with the 
equatorward migration of the magnetic fiux belts in the 
Sun from which sunspots emerge. On the other hand, at 
the equator this migration corresponds to prograde rota- 
tion which is a clear effect seen in the Sun where magnetic 
trace rs are se en to rotate faster, i.e., in the prograde direc- 
tion (|Gizon et al,, 20 03). Eve n sunspots rot ate faster than 
the gas itself (jPulkkinen fc Tuominenl . I1998D . 

One of our goals for future work is to verify the present 
findings in direct numerical simulations. Such simulations 
would also allow us to determine new turbulent transport 
coefficients, similar to the Qp parameter invoked in the 
present study. Such additional parameters yield new effects, 
some of which could be important for applications to the 
Sun. 

Finally, let us comment of the issue of scale separation. 
As discussed above, in solar mixing length theory, the cor- 
relation length of the turbulent eddies is expected to scale 
with the press ure scale h eight such that kfHp is constant 
and about 2.4 (jKemel et al .. 2012c). Theoretical considera- 
tions have shown further that the growth rate of NEMPI is 
proportional to kfHp. Since rotation is known to decrease 



the size of the turbulent eddies, i.e., to increase the value of 
kf, one might be tempted to speculate that rotation could 
even enhance the growth rate of NEMPI. However, in view 
of the present results, this seems now unlikely. 
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